Optimal Implementation Parameters of a Nonlinear Electrical Impedance Tomography Method Using the Complete Electrode Model

This study discusses a nonlinear electrical impedance tomography (EIT) technique under different analysis conditions to propose its optimal implementation parameters. The forward problem for calculating electric potential is defined by the complete electrode model. The inverse problem for reconstructing the target electrical conductivity profile is presented based on a partial-differential-equation-constrained optimization approach. The electrical conductivity profile is iteratively updated by solving the Karush–Kuhn–Tucker optimality conditions and using the conjugate gradient method with an inexact line search. Various analysis conditions such as regularization scheme, number of electrodes, current input patterns, and electrode arrangement were set differently, and the corresponding results were compared. It was found from this study that the proposed EIT method yielded appropriate inversion results with various parameter settings, and the optimal implementation parameters of the EIT method are presented. This study is expected to expand the utility and applicability of EIT for the non-destructive evaluation of structures.


Introduction
Electrical impedance tomography (EIT) is a non-destructive evaluation method through which the electrical properties of part of a structure are determined using measured data from surface electrodes. This method is highly applicable in medical imaging, industrial process monitoring, and geotechnical site characterization because of its ease of use in field experimentation, economic feasibility, and superior ability to penetrate the target. For example, the EIT method coupled with convolutional neural networks has been explored for reconstructing human organ boundaries [1], and EIT using an optimal control theory has been developed for pathological diagnoses such as cancer detection [2]. In industrial processes, EIT can be a useful tool to monitor the mixing process of chemical materials [3,4] and evaluate the dredging process's condition in real-time by monitoring the material flowing through the dredging pipe [5,6]. Recent applications of EIT in civil engineering include the characterization of layered soils [7], crack detection in pipes buried in the ground [8], and ground contamination monitoring for remediation strategies [9].
EIT has many advantages; however, it requires improvements in its mathematical modeling, numerical analysis, and implementation techniques to increase the accuracy of solutions and to broaden its application scope to further extensive fields. Previous studies on EIT suggest that the quality of the inverse tomographic images obtained using measured electric potentials is sensitive to electrode arrangement and current input patterns [10][11][12][13][14][15][16]. Graham et al. [11] presented seven electrode placement configurations (planar, planar-offset, planar-opposite, zigzag, zigzag-offset, zigzag-opposite, square) in which the electrodes were arranged in two parallel planes of eight electrodes each, with electrodes equispaced around a medium. These configurations were applied to three-dimensional (3D) EIT, and calculated electric potential and experimentally obtained data is smaller than that com puted using other models, such as a point electrode model. Figure 1 shows the configura tion of a two-dimensional (2D) square domain with electrodes on its surface. where Ω denotes structural domain, is the scalar-valued electric potential to be calcu lated, is electrical conductivity, is the outward unit normal to the boundary ∂Ω, Γ is the th electrode boundary, is the contact impedance of , is the injected curren at , is the electric potential at to be calculated, and is the number of electrodes Equation (1) is a Laplace equation for the electric potential ( ). Equation (2) is a Robin type boundary condition describing the electric potential at Γ , and Equations (3) and (4 are Neumann boundary conditions for ( ). To ensure the existence and uniqueness o the solution, the following continuity condition, = 0, (5 is added to the model. For setting the reference point of the electric potential, = 0 (6 must be satisfied. For the variational form of the boundary value problem, Equation (1) is multiplied by a test function ( ) ∈ (Ω), and then integrated over the domain Ω using boundary conditions (2), (3), and (4). Meanwhile, Equation (2) is integrated ove Γ , multiplied by a test value , and then summed for all electrodes. Adding the tw equations results in a variational form, expressed as follows [24]: (7 Introducing finite element approximations to the electric potential ( ) and the tes function ( ) results in a linear system of equations, where the electric potential a The CEM for calculating the electric potential due to the current input can be expressed as a boundary value problem, as follows: u + z l σ ∂u ∂n = U l , l = 1, 2, · · · , L, Γ E l σ ∂u ∂n dΓ = I l l = 1, 2, · · · , L, where Ω denotes structural domain, u is the scalar-valued electric potential to be calculated, σ is electrical conductivity, n is the outward unit normal to the boundary ∂Ω, Γ E l is the lth electrode boundary, z l is the contact impedance of E l , I l is the injected current at E l , U l is the electric potential at E l to be calculated, and L is the number of electrodes. Equation (1) is a Laplace equation for the electric potential u(x). Equation (2) is a Robin-type boundary condition describing the electric potential at Γ E l , and Equations (3) and (4) are Neumann boundary conditions for u(x). To ensure the existence and uniqueness of the solution, the following continuity condition, is added to the model. For setting the reference point of the electric potential, must be satisfied. For the variational form of the boundary value problem, Equation (1) is multiplied by a test function v(x) ∈ H 1 (Ω), and then integrated over the domain Ω using boundary conditions (2), (3), and (4). Meanwhile, Equation (2) is integrated over Γ E l , multiplied by a test value V l , and then summed for all electrodes. Adding the two equations results in a variational form, expressed as follows [24]: Introducing finite element approximations to the electric potential u(x) and the test function v(x). results in a linear system of equations, where the electric potential u at Sensors 2022, 22, 6667 4 of 28 each node and the electrode potential U l can be calculated. The stiffness matrix and the right-hand-side vector of the linear system can be found in [19,20].

Setting for Numerical Analysis
The forward CEM solution is validated by comparing it to the solution obtained by ANSYS Mechanical APDL. Figure 2a shows the configuration of a homogeneous square domain with a side length of 10 m. The electrical conductivity of the domain is 0.01 S/cm. A total of 10 electrodes are distributed equally on the left and right sides, as shown in the figure. Figure 2b shows a finite element model consisting of 1600 eight-node square elements with sides of 0.25 m. The length of one electrode is 0.25 m, equal to the size of the finite element. A current of 0.1 A is input into the electrodes placed on the left side of the domain, and it flows out through the electrodes placed on the right side. The contact impedance of the electrodes is 1 × 10 −5 Ω·m 2 . each node and the electrode potential can be calculated. The stiffness matrix and the right-hand-side vector of the linear system can be found in [19,20].

Setting for Numerical Analysis
The forward CEM solution is validated by comparing it to the solution obtained by ANSYS Mechanical APDL. Figure 2a shows the configuration of a homogeneous square domain with a side length of 10 m. The electrical conductivity of the domain is 0.01 S/cm. A total of 10 electrodes are distributed equally on the left and right sides, as shown in the figure. Figure 2b shows a finite element model consisting of 1600 eight-node square elements with sides of 0.25 m. The length of one electrode is 0.25 m, equal to the size of the finite element. A current of 0.1 A is input into the electrodes placed on the left side of the domain, and it flows out through the electrodes placed on the right side. The contact impedance of the electrodes is 1 × 10 Ω • m .  each node and the electrode potential can be calculated. The stiffness matrix and the right-hand-side vector of the linear system can be found in [19,20].

Setting for Numerical Analysis
The forward CEM solution is validated by comparing it to the solution obtained by ANSYS Mechanical APDL. Figure 2a shows the configuration of a homogeneous square domain with a side length of 10 m. The electrical conductivity of the domain is 0.01 S/cm. A total of 10 electrodes are distributed equally on the left and right sides, as shown in the figure. Figure 2b shows a finite element model consisting of 1600 eight-node square elements with sides of 0.25 m. The length of one electrode is 0.25 m, equal to the size of the finite element. A current of 0.1 A is input into the electrodes placed on the left side of the domain, and it flows out through the electrodes placed on the right side. The contact impedance of the electrodes is 1 × 10 Ω • m .

PDE-Constrained Optimization
The inverse medium problem for reconstructing the electrical conductivity profile of a structure using measured electric potential data at surface electrodes can be presented as the following PDE-constrained optimization problem: The objective functional comprises a misfit functional and a regularization term. The misfit functional is expressed as the sum of the squared differences of the calculated electric potential and the measured electric potential at electrode . This optimization problem is constrained by Equation (1), which is the governing equation of the CEM, and boundary conditions (2)-(4). To relieve the ill-posedness present in such an inverse problem, the regularization term ( ) for the electrical conductivity is included in the objective functional .

Regularization Schemes
In this study, Tikhonov (TN) [25] and total variation (TV) [26] regularization schemes are used to investigate the regularization effect. For TN regularization, the regularization term ( ) can be expressed by Equation (9):

PDE-Constrained Optimization
The inverse medium problem for reconstructing the electrical conductivity profile of a structure using measured electric potential data at surface electrodes can be presented as the following PDE-constrained optimization problem: The objective functional J comprises a misfit functional and a regularization term. The misfit functional is expressed as the sum of the squared differences of the calculated electric potential U l and the measured electric potential U m l at electrode E l . This optimization problem is constrained by Equation (1), which is the governing equation of the CEM, and boundary conditions (2)-(4). To relieve the ill-posedness present in such an inverse problem, the regularization term γ(σ) for the electrical conductivity σ is included in the objective functional J.

Regularization Schemes
In this study, Tikhonov (TN) [25] and total variation (TV) [26] regularization schemes are used to investigate the regularization effect. For TN regularization, the regularization term γ(σ) can be expressed by Equation (9): For TV regularization, γ(σ) can be expressed by Equation (10): where R σ is a regularization factor that controls the penalty for the spatial variation of electrical conductivity σ(x). In Equation (10), a small parameter β is included to make γ TV (σ) differentiable when ∇σ = 0. Generally, it is expected that TN regularization would be suitable for reconstructing a smooth target profile. On the other hand, TV regularization is expected to perform better when reconstructing a sharply varying target profile.

First-Order Optimality Conditions
The Lagrange multiplier method was used to convert the PDE-constrained optimization problem written in Equation (8) into an unconstrained optimization problem. The objective functional J can be augmented using Equations (1) and (2) to construct the Lagrangian functional L: where w and W l are Lagrange multipliers multiplied to the left-hand terms of the governing equation and boundary conditions, respectively. The electrical conductivity σ(x) that minimizes the Lagrangian L is the solution to the inverse problem. For the optimal solution to this problem, the first-order optimality conditions of the Lagrangian are enforced. In other words, the first variation of L with respect to adjoint variables w and W l , state variables u and U l , and control variable σ is enforced to vanish. There results the state problem for u and U l , the adjoint problem for w and W l , and the control problem for σ, respectively. Solving the three problems simultaneously in the reduced space of the control variable yields the optimal solution of the material profile σ(x).

First Optimality Condition: State Problem
The state equation and the corresponding boundary conditions can be obtained from the stationarity requirement that the first variation of the Lagrangian with respect to adjoint variables w and W l must be 0 (δ w L = 0, δ W l L = 0, l = 1, 2, · · · L). The derived state problem is identical to the forward problem in Equations (1)-(4).

Second Optimality Condition: Adjoint Problem
The adjoint equation and the corresponding boundary conditions can be obtained from the stationarity requirement that the first variation of the Lagrangian with respect to state variables u and U l must be 0 (δ u L = 0, δ U l L = 0, l = 1, 2, · · · L). The derived adjoint problem can be described as follows [19,20]: w + z l σ ∂w ∂n = −W l , l = 1, 2, · · · , L, x ∈ Γ E l , Equation (12) is the governing equation for the adjoint variable w(x). It has a differential operator similar to (1), the state equation. Equations (13)- (15) are the boundary conditions of the adjoint problem. Equation (14) indicates the source of the adjoint problem, which depends on the misfit of the electric potential at electrodes. The adjoint problem can also be solved by the finite element method in a manner similar to the state problem [20].

Third Optimality Condition: Control Problem
The control equation and the corresponding boundary conditions can be obtained from the stationarity requirement that the first variation of the Lagrangian with respect to the control variable σ must be 0 (δ σ L = 0). The derived control problem can be described as follows [19,20]: In deriving Equation (16), the TN regularization scheme was used. If the TV scheme were used instead, Equation (16) could be replaced by where H is the Hessian matrix of σ(x). It is different from the Hessian for the Gauss-Newton inversion, which consists of the second Fréchet derivatives of the Lagrangian [27][28][29]. The solution σ(x) of the control problem can be calculated once the state and adjoint solutions u, w, and W l are obtained. The state, adjoint, and control problems derived from the firstorder optimality conditions of the Lagrangian indicate the Karush-Kuhn-Tucker (KKT) conditions for this optimization problem.

Material Property Update
The first and second optimality conditions are satisfied by solving the state and adjoint problems, respectively. Because only the true profile of σ(x) exactly satisfies the control problem, the material profile σ(x) must be updated to satisfy the third optimality condition. The procedure of updating the control variable σ(x) using the state and adjoint solutions is as follows:

1.
Assume the initial electrical conductivity profile of a structure to be investigated, then calculate the electric potential u and U l due to the current input through the surface electrodes.

2.
Calculate the adjoint solutions w and W l using the state solution U l . 3.
Using the state and adjoint solutions, calculate the gradient of the Lagrangian with respect to the control variable σ, as follows: In Equation (20), the TN regularization scheme was used. If the TV scheme were assumed, then the Lagrangian gradient for σ would be

4.
Update the electrical conductivity at each node using a line search method. Equations (17) and (18) are not precisely enforced in updating the electrical conductivity at boundaries since they are complicated to implement. Instead, one can enforce that the normal derivative of σ(x) be zero along the boundary for computational simplicity.

Conjugate Gradient Method with an Inexact Line Search
The search direction for the optimal solution of the control variable σ is determined using the Flétcher-Reeves conjugate gradient method with an inexact line search. Let g k denote the discrete reduced gradient at the kth inversion iteration.
Thereafter, the electrical conductivity vector σ k comprising nodal values of σ is updated via where d k is the search direction vector at σ k , and α. is the step length in the direction of d k . The step length can be determined by a backtracking algorithm in Table 1 [30]. In this work, ρ = 0.5 was used. Table 1. Backtracking algorithm to determine the step length α.

Regularization Factor Continuation Scheme
The choice of regularization factor R σ in Equations (20) and (21) considerably affects the reconstruction of the electrical conductivity profile because it controls the amount of imposed penalty on high-frequency oscillations of the material properties. In this study, a regularization factor continuation scheme [30][31][32] was used to determine the optimal regularization factor at each inversion iteration. The reduced gradients in Equations (20) and (21) can be rewritten as where R σ ∇ σ J r denotes the gradient of the regularization functional and ∇ σ J m the gradient of the misfit functional. In the case of the TN regularization, If the TV regularization scheme were used, The first term of Equation (24), R σ (∇ σ J r ), penalizes spatial oscillations in the reconstructed profile, such that a higher R σ results in a smoother reconstructed profile. A balance between these two terms can be imposed using Therefore, R σ can be calculated at each iteration as where ε is a weight factor, which plays the role of regularization effect controller. ε = 1 results in the maximum regularization effect, and ε = 0 indicates no regularization.

Numerical Studies for Optimal EIT Parameters
Consider a square domain with a side length of 10 m, as shown in Figure 1, surrounded by electrodes with a contact impedance of 1 × 10 −5 Ω·m 2 . The initial value of the regularization factor R σ is 1.0, and the weight factor ε for the regularization factor continuation scheme is set to 0.5. The parameter β for TV regularization is assumed to be This study compares the inversion results according to the change in various implementation parameters using a response misfit and a relative L 2 -error. The response misfit F m , as part of the objective functional J in Equation (8), can be written as The relative L 2 -error of electrical conductivity, E L 2 can be written as where A is the total area of the domain, σ tg is the target electrical conductivity, and σ inv is the reconstructed electrical conductivity.

Regularization Effect
For evaluating the regularization effect on the EIT, the TN and TV regularization schemes are explored in the inversion for a three-layer heterogeneous medium. Figure 5 shows the target electrical conductivity profile with three layers and an initial guess for inversion. The target values of the electrical conductivity are 0.01 S/m for y ≥ −4 m, 0.03 S/m for −7 m ≤ y < −4 m, and 0.05 S/m for y < −7 m. The values are typical of air-dried concrete materials in various conditions [33]. Examples of the profile heterogeneity in Figure 5 include fiber-reinforced composite sandwich plates and concrete specimens under curing. The initial guess for the inversion σ ini = 0.03 S/m. 40 electrodes are arranged on all sides of the square medium with equal spacing; thereafter, the current is injected into the electrodes attached on the top and left sides of the structure, and then flows out through the electrodes on the right and bottom sides. The magnitude of the current is uniform at 0.1 A. Figure 6 shows the inversion results of the three-layer electrical conductivity profile at 3000 iterations using TN and TV regularization schemes. When the TV scheme was used, the target profile was reconstructed clearly and stably, especially at the interface of the layers. This shows that the TV scheme performs well in the EIT framework when reconstructing a sharply varying profile. Figure 7 shows the response misfit and the relative L 2 -error E L 2 against iteration numbers. The misfit is reduced by 99.8% from its initial value for the TN scheme, and by 99.9% for the TV at 1500 iterations. Figure 8 shows the measured, initial, and calculated electric potentials in the inversion. After the inversion, the calculated potential values nearly coincide with the measured values, indicating the successful reconstruction of the target profile. Figure 9 shows the target and reconstructed conductivity profiles at y = −7.  As mentioned in Section 3.6, the regularization factor considerably affects the reconstruction of the electrical conductivity profile. Figure 10 shows the change in the regularization factor during the inversion using the regularization factor continuation  As mentioned in Section 3.6, the regularization factor considerably affects the reconstruction of the electrical conductivity profile. Figure 10 shows the change in the regularization factor during the inversion using the regularization factor continuation As mentioned in Section 3.6, the regularization factor R σ considerably affects the reconstruction of the electrical conductivity profile. Figure 10 shows the change in the regularization factor during the inversion using the regularization factor continuation scheme. The value of the weight factor ε was assumed to be 0.5, 0.3, 0.1, and 0.05. The higher the weight factor ε, the larger the R σ value. It can also be seen that the value of R σ fluctuates significantly in the latter part of the inversion as the inverted profile approaches the target. scheme. The value of the weight factor was assumed to be 0.5, 0.3, 0.1, and 0.05. The higher the weight factor , the larger the value. It can also be seen that the value of fluctuates significantly in the latter part of the inversion as the inverted profile approaches the target. For demonstrating the effectiveness of the regularization factor continuation scheme, the inversion tried fixed regularization factors in the same setting. Figures 11 and 12 show the reconstructed three-layer electrical conductivity profiles using different fixed regularization factors. The layer interfaces are not properly recovered when is large, but the stratum is better reconstructed when the factor is small. Again, the TV scheme yielded sharper profile reconstruction than the TN. Compared to using fixed regularization factors, the continuation scheme determines the regularization factor adaptively at each inversion iteration, making it possible to skip multiple inversion attempts to find the optimal fixed regularization factor. Figure 13 shows the variation of the relative -error in Equation (32) to iteration numbers. In the case of the continuous regularization factor, the relative -error is similar to fixed regularization factor cases at the early inversion stage, but eventually becomes smaller as the inversion progresses.  For demonstrating the effectiveness of the regularization factor continuation scheme, the inversion tried fixed regularization factors in the same setting. Figures 11 and 12 show the reconstructed three-layer electrical conductivity profiles using different fixed regularization factors. The layer interfaces are not properly recovered when R σ is large, but the stratum is better reconstructed when the factor is small. Again, the TV scheme yielded sharper profile reconstruction than the TN. Compared to using fixed regularization factors, the continuation scheme determines the regularization factor adaptively at each inversion iteration, making it possible to skip multiple inversion attempts to find the optimal fixed regularization factor. Figure 13 shows the variation of the relative L 2 -error in Equation (32) to iteration numbers. In the case of the continuous regularization factor, the relative L 2 -error is similar to fixed regularization factor cases at the early inversion stage, but eventually becomes smaller as the inversion progresses. scheme. The value of the weight factor was assumed to be 0.5, 0.3, 0.1, and 0.05. The higher the weight factor , the larger the value. It can also be seen that the value of fluctuates significantly in the latter part of the inversion as the inverted profile approaches the target. For demonstrating the effectiveness of the regularization factor continuation scheme, the inversion tried fixed regularization factors in the same setting. Figures 11 and 12 show the reconstructed three-layer electrical conductivity profiles using different fixed regularization factors. The layer interfaces are not properly recovered when is large, but the stratum is better reconstructed when the factor is small. Again, the TV scheme yielded sharper profile reconstruction than the TN. Compared to using fixed regularization factors, the continuation scheme determines the regularization factor adaptively at each inversion iteration, making it possible to skip multiple inversion attempts to find the optimal fixed regularization factor. Figure 13 shows the variation of the relative -error in Equation (32) to iteration numbers. In the case of the continuous regularization factor, the relative -error is similar to fixed regularization factor cases at the early inversion stage, but eventually becomes smaller as the inversion progresses.

Parametric Studies for Optimal EIT Result
In this study, three different analysis conditions were explored to derive optimal parameters for performing the EIT in heterogeneous media. The conditions considered are the number of electrodes, spatial current input pattern, and electrode arrangement. Figure  14 shows the target electrical conductivity profile with two layers and the initial guess for inversion. The target values of the electrical conductivity are 0.02 S/m for ≥ −5 m and 0.04 S/m for < −5 m. The initial guess of the profile is homogeneous, with = 0.03 S/m. The TV regularization with the regularization factor continuation scheme is used in all parametric studies.

Parametric Studies for Optimal EIT Result
In this study, three different analysis conditions were explored to derive optimal parameters for performing the EIT in heterogeneous media. The conditions considered are the number of electrodes, spatial current input pattern, and electrode arrangement. Figure  14 shows the target electrical conductivity profile with two layers and the initial guess for inversion. The target values of the electrical conductivity are 0.02 S/m for ≥ −5 m and 0.04 S/m for < −5 m. The initial guess of the profile is homogeneous, with = 0.03 S/m. The TV regularization with the regularization factor continuation scheme is used in all parametric studies.

Parametric Studies for Optimal EIT Result
In this study, three different analysis conditions were explored to derive optimal parameters for performing the EIT in heterogeneous media. The conditions considered are the number of electrodes, spatial current input pattern, and electrode arrangement. Figure 14 shows the target electrical conductivity profile with two layers and the initial guess for inversion. The target values of the electrical conductivity are 0.02 S/m for y ≥ −5 m and 0.04 S/m for y < −5 m. The initial guess of the profile is homogeneous, with σ ini = 0.03 S/m. The TV regularization with the regularization factor continuation scheme is used in all parametric studies.

Parametric Studies for Optimal EIT Result
In this study, three different analysis conditions were explored to derive optimal parameters for performing the EIT in heterogeneous media. The conditions considered are the number of electrodes, spatial current input pattern, and electrode arrangement. Figure  14 shows the target electrical conductivity profile with two layers and the initial guess for inversion. The target values of the electrical conductivity are 0.02 S/m for ≥ −5 m and 0.04 S/m for < −5 m. The initial guess of the profile is homogeneous, with = 0.03 S/m. The TV regularization with the regularization factor continuation scheme is used in all parametric studies.

Number of Electrodes
The number of electrodes is expected to significantly affect the inversion result because it impacts the electric potential field and the amount of measured potential data. In this work, 8,20,40, and 80 electrodes are tested for inversion, as shown in Figure 15. An equal number of electrodes are placed on each side of the square medium. The current is supplied into the electrodes attached to the top and left sides of the medium, and then flows out from the electrodes on the right and bottom sides. The magnitude of the current is uniform at 0.1 A.

Number of Electrodes
The number of electrodes is expected to significantly affect the inversion result because it impacts the electric potential field and the amount of measured potential data. In this work, 8,20,40, and 80 electrodes are tested for inversion, as shown in Figure 15. An equal number of electrodes are placed on each side of the square medium. The current is supplied into the electrodes attached to the top and left sides of the medium, and then flows out from the electrodes on the right and bottom sides. The magnitude of the current is uniform at 0.1 A.  Figure 16 shows the reconstructed two-layer electrical conductivity profiles at 2000 inversion iterations using the presented number of electrodes. In the case of eight electrodes, the target profile is not reconstructed properly compared to other cases because the number of electrodes on the surface is significantly insufficient. In other cases, the target profile is reasonably reconstructed. Figure 17a shows the variation of response misfit to iteration numbers during the inversion. In the case of eight electrodes, the misfit ℱ changes extremely unstable, even for small iteration numbers. In other cases, ℱ decreases by a factor of 10 to 10 . Figure 17b exhibits the variation of the relative -  Figure 16 shows the reconstructed two-layer electrical conductivity profiles at 2000 inversion iterations using the presented number of electrodes. In the case of eight electrodes, the target profile is not reconstructed properly compared to other cases because the number of electrodes on the surface is significantly insufficient. In other cases, the target profile is reasonably reconstructed. Figure 17a shows the variation of response misfit to iteration numbers during the inversion. In the case of eight electrodes, the misfit F m changes extremely unstable, even for small iteration numbers. In other cases, F m decreases by a factor of 10 −3 to 10 −4 . Figure 17b exhibits the variation of the relative L 2 -error to iteration numbers. In the case of 20 electrodes, the error is smaller than 40 and 80 electrode cases after about 750 iterations. Thus, a higher number of electrodes does not necessarily increase the inversion quality. Figure 18 shows the measured, initial, and calculated electric potentials in the inversion using four different electrode numbers. After the inversion, the calculated potential values coincide with the measured values. Figure 19 shows the reconstructed electrical conductivity profiles at y = −7.5 m, y = −2.5 m, x = −3.0 m, and x = 3.0 m after the inversion using four different electrode numbers. Except for the case of eight electrodes, the calculated conductivity values capture the target sufficiently well in all locations. Figure 20 presents the inverted profiles using TN and TV regularization schemes in the case of eight surface electrodes. In comparison with the TV scheme, the inversion result has been improved in the case of TN regularization. Therefore, it is more appropriate to use TN regularization when the number of electrodes is small. error to iteration numbers. In the case of 20 electrodes, the error is smaller than 40 and 80 electrode cases after about 750 iterations. Thus, a higher number of electrodes does not necessarily increase the inversion quality. Figure 18 shows the measured, initial, and calculated electric potentials in the inversion using four different electrode numbers. After the inversion, the calculated potential values coincide with the measured values. Figure  19 shows the reconstructed electrical conductivity profiles at = −7.5 m, = −2.5 m, = −3.0 m, and = 3.0 m after the inversion using four different electrode numbers. Except for the case of eight electrodes, the calculated conductivity values capture the target sufficiently well in all locations. Figure 20 presents the inverted profiles using TN and TV regularization schemes in the case of eight surface electrodes. In comparison with the TV scheme, the inversion result has been improved in the case of TN regularization. Therefore, it is more appropriate to use TN regularization when the number of electrodes is small.

Current Input Pattern
Uniform and cosine input patterns are investigated as the current input pattern for inversion. The uniform pattern is the same as that used in Section 4.2.1. In the case of the cosine pattern, this study introduced four current input phases to the inversion. Equations (33) and (34) describe the uniform and cosine current input patterns, respectively.   In Equation (34), N denotes the number of electrodes, and k is a specific electrode number. In this numerical experiment, N is 40, as shown in Figure 15c. Figure 21 shows the inverted two-layer electrical conductivity profiles at 2000 inversion iterations using the described current input patterns. Despite some differences in the results, especially at the layer interface, all current patterns successfully reconstructed the target profile. Figure 22 shows the response misfit and relative L 2 -error to iteration numbers in the inversion using the current input patterns. In all cases, the misfit is reduced by more than 99.9% from the initial misfit at 500 iterations. Figure 22b shows that the relative L 2 -error in the case of the uniform pattern is larger than that of the cosine patterns. The error is smallest when the phase α = 0. Figure 23 shows the measured, initial, and calculated electric potentials in the inversion using the current input patterns. Again, the calculated potential values are almost identical to the measured values. Figure 24 shows the reconstructed conductivity profiles at y = −7.5 m, y = −2.5 m, x = −3.0 m, and x = 3.0 m after the inversion using the current input patterns. All patterns reconstruct the target electrical conductivity profile fairly well.
In Equation (34), denotes the number of electrodes, and is a specific electrode number. In this numerical experiment, is 40, as shown in Figure 15c. Figure 21 shows the inverted two-layer electrical conductivity profiles at 2000 inversion iterations using the described current input patterns. Despite some differences in the results, especially at the layer interface, all current patterns successfully reconstructed the target profile. Figure  22 shows the response misfit and relative -error to iteration numbers in the inversion using the current input patterns. In all cases, the misfit is reduced by more than 99.9% from the initial misfit at 500 iterations. Figure 22b shows that the relative -error in the case of the uniform pattern is larger than that of the cosine patterns. The error is smallest when the phase = 0. Figure 23 shows the measured, initial, and calculated electric potentials in the inversion using the current input patterns. Again, the calculated potential values are almost identical to the measured values. Figure 24 shows the reconstructed conductivity profiles at = −7.5 m, = −2.5 m, = −3.0 m, and = 3.0 m after the inversion using the current input patterns. All patterns reconstruct the target electrical conductivity profile fairly well.      Figure 25 shows two different types of electrode arrangements used for the EIT. The first arrangement type is to place electrodes on all sides of the medium, and the second is to position them only on two sides. The total number of electrodes is 40. The uniform and cosine (α = π 2 ) current input patterns are used for this case. Figure 26 shows the inversion results at 2000 iterations using the two electrode arrangements. The target profile has been reconstructed well for both current input patterns. However, the quality of reconstruction at the layer interface is better for the all-side arrangement. This happens when the electrodes are attached to only two sides of the structure, and the inner information of the top and bottom parts cannot be sufficiently captured by surface electrodes. Figure 27 shows the response misfit and the relative L 2 -error to iteration numbers during the inversion using the two electrode arrangements. For all cases, the misfit decreased by a factor of 10 −3 to 10 −4 . As shown in Figure 27b, the relative L 2 -error for the all-side arrangement is smaller than for the two-side arrangement. In addition, the error for the cosine current pattern is smaller than for the uniform pattern. Figure 28 shows the measured, initial, and calculated electric potentials in the inversion using the two electrode arrangements. The excellent agreement of the calculated and measured electric potential values demonstrates the feasibility of the inversion. Figure 29 presents the target and reconstructed electrical conductivity profiles at y = −7.5 m, y = −2.5 m, x = −3.0 m, and x = 3.0 m. Overall, the all-side arrangement results in better reconstruction of the target profile than the two-side arrangement. Figure 25 shows two different types of electrode arrangements used for the EIT. The first arrangement type is to place electrodes on all sides of the medium, and the second is to position them only on two sides. The total number of electrodes is 40. The uniform and cosine ( = ) current input patterns are used for this case.  Figure 26 shows the inversion results at 2000 iterations using the two electrode arrangements. The target profile has been reconstructed well for both current input patterns. However, the quality of reconstruction at the layer interface is better for the all-side arrangement. This happens when the electrodes are attached to only two sides of the structure, and the inner information of the top and bottom parts cannot be sufficiently captured by surface electrodes. Figure 27 shows the response misfit and the relative -error to iteration numbers during the inversion using the two electrode arrangements. For all cases, the misfit decreased by a factor of 10 to 10 . As shown in Figure 27b, the relativeerror for the all-side arrangement is smaller than for the two-side arrangement. In addition, the error for the cosine current pattern is smaller than for the uniform pattern. Figure  28 shows the measured, initial, and calculated electric potentials in the inversion using the two electrode arrangements. The excellent agreement of the calculated and measured electric potential values demonstrates the feasibility of the inversion. Figure      (c) (d)  Figure 30 shows the variation of ‖ ‖ corresponding to the number of electrodes in the inversion using the all-side electrode arrangement. The values of ‖ ‖ were calculated after 2000 inversion iterations. As the number of electrodes increases, the relative -error tends to decrease. However, when the number of electrodes is 40 or more, there is a slight difference in the error. In addition, the error reduces when the cosine current input pattern is used, especially when = 0 or = /2. A similar error trend can be observed in the case of the two-side electrode arrangement. Figure 31 presents the variation of ‖ ‖ corresponding to the number of electrodes in the inversion cases of different electrode arrangement and current input pattern. The relative -error is smaller in the case of the all-side electrode arrangement than in the two-side arrangement. Tables 2 and 3 show the relative -error and the relative misfit (|ℱ | /|ℱ | ) for all cases of the number of electrodes, current input pattern, and electrode arrangement discussed so far. The misfit, |ℱ | , is the one immediately before the start of the misfit oscillation, and |ℱ | is the initial misfit. From the error values in the tables, one can choose the optimal implementation parameters of the described EIT. The first parameter set for which the value of ‖ ‖ is minimal is 80 electrodes, cosine current input pattern with = /2, and the all-side electrode arrangement. The second parameter set for which the value of |ℱ | /|ℱ | is minimal is the same as the first set except for the number of electrodes, which is 40. The minimum error values are shaded in the tables.  Figure 30 shows the variation of E L 2 corresponding to the number of electrodes in the inversion using the all-side electrode arrangement. The values of E L 2 were calculated after 2000 inversion iterations. As the number of electrodes increases, the relative L 2 -error tends to decrease. However, when the number of electrodes is 40 or more, there is a slight difference in the error. In addition, the error reduces when the cosine current input pattern is used, especially when α = 0 or α = π/2. A similar error trend can be observed in the case of the two-side electrode arrangement. Figure 31 presents the variation of E L 2 corresponding to the number of electrodes in the inversion cases of different electrode arrangement and current input pattern. The relative L 2 -error is smaller in the case of the all-side electrode arrangement than in the two-side arrangement. Tables 2 and 3 show the relative L 2 -error and the relative misfit (|F m | opt /|F m | ini ) for all cases of the number of electrodes, current input pattern, and electrode arrangement discussed so far. The misfit, |F m | opt , is the one immediately before the start of the misfit oscillation, and |F m | ini is the initial misfit. From the error values in the tables, one can choose the optimal implementation parameters of the described EIT. The first parameter set for which the value of E L 2 is minimal is 80 electrodes, cosine current input pattern with α = π/2, and the all-side electrode arrangement. The second parameter set for which the value of |F m | opt /|F m | ini is minimal is the same as the first set except for the number of electrodes, which is 40. The minimum error values are shaded in the tables. oscillation, and |ℱ | is the initial misfit. From the error values in the tables, one can choose the optimal implementation parameters of the described EIT. The first parameter set for which the value of ‖ ‖ is minimal is 80 electrodes, cosine current input pattern with = /2, and the all-side electrode arrangement. The second parameter set for which the value of |ℱ | /|ℱ | is minimal is the same as the first set except for the number of electrodes, which is 40. The minimum error values are shaded in the tables.        Figure 32 shows the reconstructed three-layer profiles using the optimal implementation parameters. The target and initial guess of the electrical conductivity profile are the same as those in Figure 5. The inverted profiles are obtained at 3000 iterations. The quality of profile reconstruction using the optimal parameter set is considerably better than the result shown in Figure 6, especially at the layer interface. Figure 33 shows the misfit and the relative L 2 -error against iteration numbers in the inversion using the optimal parameter sets. The relative L 2 -error decreased by 95.1% compared to the initial value when the first parameter set is used. It also decreased by 94.4% when using the second parameter set. The reduction rate of E L 2 is slightly lower than in Figure 7b. first parameter set is used. It also decreased by 94.4% when using the second parameter set. The reduction rate of ‖ ‖ is slightly lower than in Figure 7b.

Conclusions
This study investigated optimal implementation parameters for a nonlinear EIT technique using the CEM. The EIT method is based on PDE-constrained optimization, which reconstructs the electrical conductivity profile by solving the KKT conditions iteratively. first parameter set is used. It also decreased by 94.4% when using the second parameter set. The reduction rate of ‖ ‖ is slightly lower than in Figure 7b.

Conclusions
This study investigated optimal implementation parameters for a nonlinear EIT technique using the CEM. The EIT method is based on PDE-constrained optimization, which reconstructs the electrical conductivity profile by solving the KKT conditions iteratively. By applying various analysis conditions, the optimal set of parameters that minimize rel-

Conclusions
This study investigated optimal implementation parameters for a nonlinear EIT technique using the CEM. The EIT method is based on PDE-constrained optimization, which reconstructs the electrical conductivity profile by solving the KKT conditions iteratively. By applying various analysis conditions, the optimal set of parameters that minimize relative L 2 -error or relative misfit in the EIT has been derived. The quality of the reconstructed profile using the optimal implementation parameters is superior to the results using a conventional parameter set.

1.
The layered profile was reconstructed more clearly when using the TV regularization scheme than TN, especially at the interface of layers. The inversion result was improved when using the regularization factor continuation scheme rather than the fixed method.

2.
A higher number of electrodes did not necessarily improve the inversion results. In addition, the TN regularization scheme produced relevant results when the number of electrodes was small.

3.
The layered profiles were successfully reconstructed for all the presented current patterns. The relative L 2 -error was smaller when the cosine pattern was used, especially when the phase α = 0 or α = π 2 . 4.
In the case of arranging electrodes on all sides of the square domain, the inversion result was improved compared to the case of arranging them only on two sides. 5.
The relative L 2 -error and the relative misfit are proper criteria for optimal implementation parameters. The relative L 2 -error was decreased by 95.1% from the initial value when using the first set of optimal parameters. It was also reduced by 94.4% when using the second set. The presented optimal parameter sets worked successfully in reconstructing layered electrical conductivity profiles.
This study is expected to expand the applicability of the nonlinear EIT method for the non-destructive evaluation of civil structures such as damage inspection, strength inspection of concrete under curing, fiber content inspection of fiber-reinforced concrete, etc.